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Abstract — We present equivalent circuits that model the in- 
teraction of microwave resonators and quantum systems. The 
circuit models are derived from a general interaction Hamilto- 
nian. Quantitative agreement between the simulated resonator 
transmission frequency, qubit Lamb shift and experimental data 
will be shown. We demonstrate that simple circuit models, using 
only linear passive elements, can be very useful in understanding 
systems where a small quantum system is coupled to a classical 
microwave apparatus. 

Index Terms — cavity quantum electrodynamics, circuit model- 
ing, coupled quantum-classical systems, quantum non demolition 
measurement 



I. Introduction 

CAVITY quantum electrodynamics has long been an active 
field to study the interaction of electromagnetic radiation 
with matter, which is a fundamentally important topic in 
physics _Q, 0, 0. Recently, cavity quantum electrodynamic 
experiments were performed in the microwave regime 0, 
151 , 151 , 171 , where the cavity consisted of a high-quality 
factor (Q) microwave coplanar resonator. The quantum sys- 
tem was a charge qubit built from two Josephson junctions. 
Similar experiments were proposed to bring other nanosystems 
(molecules 0, nanomagnets 0, fl__, flux qubits iflOl ifTTI 
lfl2l ) into interaction with coplanar waveguides. One of the 
most promising aspects of these experiments is that they can 
provide circuits which integrate quantum mechanical behavior 
and 'conventional' (high-frequency) components on a chip. 

It is well established how to construct circuit models from 
electromagnetic models for commonly used active or passive 
components, which obey 'classical' circuit theory. It has been 
also demonstrated how one can build equivalent circuit models 
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of basic two-state quantum systems |[T3l . The goal of the 
present paper is to show a systematic approach for building 
circuit models for an interacting resonator-quantum system, 
using linear approximations. This can be useful for developing 
quantum systems coupled via a resonator |[T4l or modeling the 
large external circuitry coupled to the quantum system. The 
possibility of this modeling approach was mentioned in [5|; to 
our knowledge, though, our paper is the first to analyze and 
exploit this model in detail. We investigate an experimental 
setup, similar to the one in lfT51 . composed of a high-Q, 
superconducting microwave resonator, electrically coupled to 
a quantum circuit (often referred as transmon-qubit [16|). 

During the past years the field of circuit quantum electro- 
dynamics went through a considerable growth. Many quan- 
tum phenomena known from optical and three dimensional 
microwave systems were also observed in circuit quantum 
electrodynamics and simulated, including systems of entan- 
gled quantum circuits IfTTI and resonators coupled via flux 
qubits ifTHl . Such quantum electrodynamic (QED) systems 
need further optimization of both the quantum bit and the 
resonator, such as increasing qubit decay times and optimizing 
resonator losses |fT9l . 

The natural application area of circuit models is the un- 
derstanding of the interaction between a relatively complex 
'classical' circuitry and a relatively simple quantum circuit. 
These models do not show new physics, but rather facilitate 
the engineering of coupled microwave-quantum systems. 

II. Device geometry 

The setup of a typical circuit quantum electrodynamic 
experiment is shown in Fig. Q] A main component is a 
superconducting coplanar waveguide resonator, displayed in 
Fig. 03 (a) (longitudinal section) and (b) (cross section). A su- 
perconducting resonator can easily reach an unloaded quality 
factor in the range of Qunioaded ~ 10 6 0, ll20l . The loaded 
quality factor of the resonators is designed by using finger- 
type capacitors in the central strip. These typically determine 
the lifetime of the resonator field (which should be high) 
and also the duration of the measurement of the nanosystem 
(which should be low in order to avoid energy relaxation and 
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Fig. 1. Geometry-sketch of the (a) longitudinal section of the resonator, 
showing the quantum system-loop that is placed near the field maximum and 
(b) the cross section of the investigated coplanar waveguide resonator. The 
simplified equivalent circuit is shown in (c), representing the ideal resonator. 
Here, the capacitors C\ and C2 correspond to the finger capacitors at the two 
ends of the resonator, and the load resistors are also included. The coupling 
to the quantum system is capacitive. 



decoherence of the quantum system). The superconducting 
resonator is practically lossless and close to resonance, its port 
to port behavior can be modeled by the circuit elements in Fig. 
[T] (c), where the LC parameters are: 



L /■ 



(1) 



Here, Zq is the resonator impedance and uj r is the resonator 
frequency. The 50 SI external environment is represented by 
the load resistors R; oa( j in Fig. Q] (c) which, together with the 
finger-type capacitors C\ w C2 in the circuit model and the 
resonator frequency to r , give the loaded quality factor of the 
resonator: Qi oa ded ~ l/(Ri oa d,CiUJ r ). 

The quantum circuit is placed in the gap of the resonator 
(see Fig.Q](a)) near the electric field maximum (magnetic field 
maximum) if it is an electrically (magnetically) interacting 
quantum system (i.e., charge qubit). Its presence slightly 
changes the values of the capacitance and inductance of the 
resonator and, thus, its resonance frequency. 

The total system can be divided into a resonator part, an 
interaction part and a quantum system part. The coupling in 
the case of electrically interacting systems is capacitive. 

Throughout this paper we will formulate our models as- 
suming electrically coupled systems. We will also, however, 
mention the parameters for the magnetic case in Section III. 
Instead of modeling different geometries, the resonator will be 
considered as a lumped circuit. For more detailed investigation 
of resonator geometries, the reader is referred to iflOl , fl2| . 

m. 

III. Method 

For the modeling of the interacting resonator-quantum sys- 
tem, we use a simplified (lossless) version of the circuit model 



in Fig.Q](c), and the coupling to the environment is neglected. 
The Hamiltonian is obtained by adding the energy stored in 
the ideal (LC) resonator, the energy stored in the interaction 
of the resonator with the quantum system, and the energy of 
the quantum system: 



ideal res. energy 
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In the above equation, the different parts are illustrated by the 
boxes in Fig. [T|(c). 

The resonator energy can be defined analogous to the 
classical energy of the LC circuit in Fig. Q] (c), however, the 
voltages and currents here are operators and their expectation 
value gives the voltage and current of the LC resonant circuit. 
The resonator root-mean square current (Irms) an d voltage 
(Vrms) E)> is related to the minimum (or vacuum) field 
in the resonator. They are further discussed in Appendix [A] 

The interaction energy also contains the resonator voltage 
term (V re s)- Here, g is related to the geometric capacitance 
of the resonator-quantum system C c rq , and a% is related to the 
voltage of the quantum system. The Cj, in the QED framework, 
are referred to as Pauli matrices and are the basic observables 
of the two state quantum system. The interaction term is a 
standard term used in quantum electrodynamics, and valid for 
every electrically interacting quantum system and resonator. 

The last term in (fJJ, containing u q represents the energy 
level splitting (AE = haj q ) of the quantum system. This part 
can only be derived in terms of quantum mechanics. 

For simulating the dynamics of the coupled system, we use 
the Liouville-von Neumann equation 
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with the Hamiltonian (0. Here, A is an operator that represents 
either one of the Pauli matrices (rji) or the resonator voltage 
(Vres) or current (I res ) operator. Based on the von Neumann 
equation, we find the following coupled differential equation 
system: 
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(6) 



(7) 



(8) 



Irms Vrms 
The above equations are a direct result of the von Neumann 
equation for the dissipationless case (if Tj, T2 — > 00 and 
7 — > 0). We can observe that the first three equations (0), (O 
and © are the widely used Bloch equations IF2TI written for a 



3 



two-level quantum system; the last two are equations for the 
dynamics of normalized resonator voltage and currents, thus 
we will call Eqs. (0|i-((8]l resonator-Bloch equations (RBEs). 
Dissipation is phenomenologically described by the decay and 
decoherence constants T\ and T2 [21], and the resonator loss, 
I/7 is added similar as in ||22l . The A^ variables are the 
expectation values of the Pauli matrices and represent the 
coherence vector characterizing the state of the quantum bit. 
The derivation of the resonator-Bloch equations is shown in 
Appendix [A] 

The only approximation used here is to neglect the cor- 
relations between the qubit and field 11221 . thus writing the 
expectation value of the normalized voltage and coherence 
vectors as a product of the individual expectation values: 

(cnVres^ = (Oi) (V res ^ = V res Xi. (9) 

An equation system similar to the resonator-Bloch equations 
was introduced by Jaynes and Cummings [22|, ||2"3l , written 
for the coupling of electric fields and dipole moments. We now 
use this model to find passive circuits describing the interaction 
of the LC resonator-quantum system. 

A. Modeling the system for large detuning 

In circuit quantum electrodynamic experiments, measure- 
ments are often performed in the so-called dispersive regime, 
where the qubit-resonator detuning is much larger then the 
coupling frequency: A ^> g. In this case, off-resonant interac- 
tion between the quantum system and the resonator does not 
switch the quantum system and A3 remains constant. Thus, 
we can neglect the oscillations of the inversion term in the 
resonator-Bloch equations (A3 = 0) and perform a linear 
approximation by keeping A3 = A3. This is valid only if the 
measurement time is much smaller than the decay time T\. 
For a large decay time T\ or short measurement time window, 
the system behaves linearly, and the resonator-Bloch equations 
reduce to the four coupled linear differential equations (we will 
call them linear RBEs): Eqs. ©, ©, Q and © with A 3 = A3. 
These can be modeled by the passive circuits shown in Fig. [2] 
and represent the interaction of (a) electric and (b) magnetic 
quantum systems, with the resonator in the dispersive regime. 
Interaction between the quantum system and resonator occurs 
through their coupling capacitance in Fig. [2] (a) or mutual 
inductance in Fig. 12(b), depending on the type of interaction 
(electric or magnetic). 

In the circuit models, the decoherence time from resonator- 
Bloch equations (0]i and (0 is represented by the parallel 
resistor with R\C q « T2; the series resistor, with L q fR 2 q w T2 
and the resonator decay rate 7, is given by R r C r in case 
of electric and L r /R r in case of magnetic interactions. The 
parameters of the circuits were found as a direct result of the 
linear resonator-Bloch equations and are summarized in Table 
H for the derivation see Appendix |B] This table summarizes 
the central result of our work. 

The physical parameters of the resonator are changed from 
the values in (HJ to the values in Table U due to the presence 
of the quantum system. The circuit models presented in this 
section provide two resonance peaks in every case, which 




Fig. 2. Circuit models representing the linear quantum mechanical interaction 
of the resonator and the qubit (i.e., A3 fs A3) in case of (a) electric 
interaction (e.g., charge qubit) and (b) magnetic interactions (e.g., flux qubit). 
The coupling between the resonator and quantum system is capacitive in (a) 
and inductive in (b). The quantum system behaves as a linear LC oscillator 
that changes the frequency of the resonator. The outcoupling capacitors, 
corresponding to the resonator finger-type gaps, have very small values and 
thus can be added to the resonant circuits that model the linear interaction. 

represent the oscillations of the coupled system. The coupling 
between the quantum system and the resonator changes the 
electromagnetic response of the resonator in a way that de- 
pends on the state (A3) of the quantum system. The quantum 
state can be read out non-destructively from the measurement 
on the resonator. In physics this is often referred to as quantum 
non-demolition measurement. 

B. Modeling the system in its ground state 

Next, we will investigate the ground state of the coupled 
resonator-quantum system. We will also compare our sim- 
ulation results mainly to measurements done in the ground 
state. In case of no applied field on the input port, the 
coupled resonator-quantum system will converge to the com- 
mon ground state [Ai, A 2 , A 3 , V res /VRMS, Ues/lRMs] = 
[0,0, -1,0,0], due to the decay times T x , T 2 and I/7 in the 
resonator-Bloch equations. In the case of a small applied probe 
field on the resonator port, the system will oscillate around 
its ground state and the inversion remains unchanged (i.e., 
A3 » -1). 

Rather than implementing the drive-fields in our resonator- 
Bloch equations, we start our simulations by taking a non- 
zero mode occupation of the resonator field. For this we take 
the initial condition: [0.04; 0; -0.999; 0; 0]. This corresponds 
to fields in the resonator-quantum system with an average 
photon number of approximately 0.01. Thus, for low fields 
in the resonator (approximately 0.01 photons [5], [15]), the 
nonlinear term A2K.es/VfjMS in the resonator-Bloch equations 
will oscillate near 0, as it is a product of two very small values. 
The inversion A3 will stay approximately constant and not 
change from its ground state even when the system has no 
detuning. As the inversion can be taken constant, the linear 
approximation done in the previous subsection holds, and we 
can take A§ w —1. 
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TABLE I 

Summary of the parameters calculated from the linear resonator-Bloch equations, if the inversion is constant (A 3 ^ A§). 




Detuning [MHz] Detuning [MHz] 



Fig. 3. Transmission frequency of the resonator as a function of detuning. 
Small circles represent experimental data from [15] and the dashed lines 
represent our circuit model. The anticrossing behavior is clearly visible and 
the separation is approximately given by the vacuum Rabi frequency at 
A = 0. The excellent agreement confirms the usability of our circuit model 
to calculate cavity transmission. 



Fig. 4. Frequency (Lamb) shift of the quantum bit as a function of 
qubit detuning. Small circles represent experimental data. Continuous lines 
represent results extracted from the numerical solutions of the resonator-Bloch 
equations. Dashed lines indicate results of our passive circuit model (Fig. [2] 

(a)) 



In the next section, we will compare our derived resonator- 
Bloch equation model and circuit models to experimental data, 
and show the validity of our approximations. 

IV. Results 

We have performed numerical and circuit simulations on a 
recently measured resonator-transmon qubit device from 1 15 1. 
We took as input parameters: resonator frequency uj r / (2tt) = 
6.44 GHz, impedance Zq — 50 O, vacuum Rabi frequency 
5 /(27r) = 266 MHz and cavity decay rate "//(2n) = 1.6 MHz. 
The dephasing time was approximately T 2 ~ 1 /is, which is 
known to be in good agreement with experimental findings 
||T6l . Further, the linewidth of the qubit coherence vector 
elements A!,A 2 is in good agreement with the experiment 
(approximately 3 MHz) in H151 . Simulations were performed 
in LTSPICE ll24l based on the presented circuit models, using 
frequency domain AC analysis. 

The transmission of the circuit in Fig. [2] (a) was calculated 
for the above parameters, setting A§ = — 1 as the system was 
in ground state. We investigated the dependence of the result 
with respect to quantum system frequency, while keeping the 
resonator frequency fixed. This was done by varying uj q which 
was directly related to the parameters of our circuit model. 

The empty circles in Fig. [3] represent data from [15| that 
has been transformed from flux to detuning coordinates so 



that our simulated results can be compared to the measured 
cavity transmission. We can see excellent agreement of the 
transmission frequency of our circuit (dashed line) with exper- 
imental data (small circles), for positive and negative detuning, 
showing the validity of the linear approximation and the excel- 
lent description of the anticrossing of the resonator-quantum 
system by the derived circuit model. At approximately A = 0, 
the peaks are separated by the vacuum Rabi frequency g/(2n). 
The transmission graphs do not contain amplitude data, but 
we mention here, that the reduction of the peak far from the 
resonator frequency (o; r /(27r)=6.44GHz) is also observed in 
our simulation. In the case of large detuning, the peaks are 
highly separated and the qubit peak gets highly reduced, while 
the resonator frequency peak approaches the original resonator 
frequency; this shows that for large detuning the two systems 
are decoupled. 

After successfully applying our circuit model to the de- 
scription of the cavity transmission, we performed an analysis 
on the quantum bit frequency shift (Lamb shift). The cavity 
transmission was calculated as a function of the detuning; 
the Lamb shift was extracted by subtracting the bare cavity 
frequency and the detuning (A) from the simulated qubit 
frequency. This was done by solving the full resonator-Bloch 
equations ©-([Ell numerically using a basic Runge-Kutta solver 
and by employing our circuit models. The simulation results 
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Fig. 5. Theoretical demonstration of our circuit model showing the state 
dependent resonator frequency shift, that opens the possibility of a quantum 
state measurement. The quantum system decay time was 20 fis, that was 
large enough not to influence the linewidth of the excited state. The shift 
is approximately g 2 /A. The slight disagreement of the peak positions results 
from the Fourier transform on the finite simulation window of the output 
voltage signal calculated with the resonator-Bloch equations. 

were compared to experimental data from |pl"5l . for a detuning 
ranging from 100 MHz to 600 MHz. The comparison was 
done by changing the frequency of the qubit and keeping the 
resonator frequency constant. 

As shown in Fig. [4] we see an excellent agreement between 
the experimental Lamb shift (small empty circles) and the 
numerical solutions of the resonator-Bloch equations (solid 
line) as well as our linear circuit model (dashed line). Thus, 
we can affirm that also the frequency shift of the quantum 
system (Lamb shift) is well described by our circuit model. 
To this point we only did simulations of the quantum sys- 
tem in its ground state. For the numerical solution of the 
full resonator-Bloch equations, we used the initial conditions 
[Al, A2, A3, V res /VRMS,Ires/ Ir.Ms] = [0.04,0,-0.999,0,0] 
as discussed in Section Illl-BI For the circuit model, we set 
A° = -1. 

We will now investigate the case, when the quantum system 
is away from its ground state, and its frequency is far detuned 
from the resonator frequency. In this case one can perform 
a non-demolition readout of the quantum state, i.e., read 
out the inversion without changing it. This can be done if 
the detuning tends to infinity (practically is large enough 
so that the inversion changes only slightly). We then know 
from Rabi's solutions that the inversion will not change (i.e., 
the resonator is unable to switch the quantum circuit during 
measurement). For an optimal quantum state readout, however, 
the state of the quantum system should not decay during 
measurement. Based on our circuit model and resonator-Bloch 
equations presented earlier, we now estimate the value of the 
decay time T\ for an optimal readout. 

The quantum state readout measurement is performed in 
the dispersive regime, as discussed in Section IIII-AI In the 
case of small decay times T\, if the quantum system is in 
its excited state, it will quickly relax to the ground state; this 



limits the measurement time and, thus, the Fourier window 
of our numerical simulations. On the other hand, for large 
decay times, the quantum-state dependent cavity shift can 
be simulated, as shown in Fig. The simulated cavity- 
pull linewidths extracted from the Fourier transform of the 
numerical resonator-Bloch equation solutions and from the 
frequency domain simulation of our circuit model agree well. 
This shows that the broadening of the peaks is mainly due to 
decoherence T2 and cavity decay (I/7) and not the qubit decay 
T\, which was not included in the circuit model. In this limit, 
the quantum-state dependent cavity pull is resolvable, and we 
find that for this case T\ > 20 /is (or more) would be required, 
so that the two linewidths approximately agree. For smaller 
decay times, the numerical simulations of the resonator-Bloch 
equations show high asymmetry of the two peaks due to the 
non-Lorentzian behavior of the spectral line corresponding to 
the excited qubit state (A3 = 1); this is due to the decay of 
the state during simulation. In order to avoid the decay of the 
quantum state the simulation window can be reduced. For a 
too small simulation window, however, the Fourier transform 
of the signal becomes inaccurate. 

V. Conclusions 

We have derived circuit models for the understanding and 
modeling of superconducting coplanar waveguide resonators 
interacting with quantum systems. They can easily be used to 
model the cavity transmission, Lamb shift and quantum non- 
demolition measurement. We find their application straightfor- 
ward in understanding experimental data and estimating decay 
times for an optimal quantum state readout. Our models can 
be extended, e.g., by adding further circuits to represent addi- 
tional qubits. Also, embedding quantum systems in large-scale 
classical circuitry is straightforward in our circuit models. 

Appendix A 

Derivation of the resonator-Bloch equations 

For the derivation of the resonator-Bloch equations, we need 
to review a few concepts widely used and well established in 
quantum optics. The starting point is to define the resonator 
voltage, similar as in J4), J3J. The voltages and currents on 
the LC circuit are defined as: 

V R l IS V res = a + a + , (10) 

J HMS^e S = -j(a - a + ). (11) 

Here the terms a + and a represent the emission and absorption 
operators. The commutator properties of these operators are 

[a, a + ] = 1, [a, a + a] = a, [a + , a + a] = — a + , 

[a, a] = [a+,a+] = [a + a,a + a] = 0. (12) 

The brackets represent a commutator relation (also used in the 
von Neumann equation (O). 

The normalization terms Vrms an d Irms represent the 
root-mean square voltage and current [4|, [5 1, on the resonator: 
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Vrms 

I RMS 



y/hu r /(2C r ), 

\J fku r / (2L r ). 



(13) 
(14) 



If we substitute the voltage and current operators ( dlPl i and 
(fTTl ) into the general Hamiltonian shown in (f2]i and use the 
commutator operations defined in (1121 . we can get the general 
Hamiltonian 



ideal res. energy 



interaction energy 




(15) 



^3 , 



widely used in cavity quantum electrodynamics Q, ED , [25 1, 

ED- 

For each Pauli matrix <Ji,(i = 1,2,3) and for the nor- 
malized current and voltage in ( fTOb and (fTTb . we used the 
von Neumann equation (01 for finding the resonator-Bloch 
equations. 

We now show the derivation of the dynamic equation for the 
resonator normalized current (©. The other resonator-Bloch 
equations (for coherence vector and resonator normalized 
voltage) can be derived analogously. First, the normalized 
current can be written as 



I 



RMS Ut±res 



d t I r 



-jh- 



RMS res 



3^ ^RMS ( ^resHtot Hf tl res 



(16) 



Using the properties of the creation/annihilation operators 
from ([T2I 1. the above von Neumann equation simplifies to 



I RM Swires = -W r (a + <2 + ) - g<J\. 



(17) 



Using the definition in (TTOb and taking the expectation values 
of the time dependent variables in the above equation, we get 
©. 

Appendix B 

Derivation of the equivalent linear circuit 

If the inversion is constant (A3 = A3) and the system does 
not decay (i.e., T\ — > 00), we can linearize the resonator- 
Bloch equations and thus get an equivalent circuit. In this case 
© is automatically fulfilled, and the other resonator-Bloch 
equations (0, (|7]i and ([H) simplify to the linear differential 
equation system: 



Ai 
A 2 

*RMS* res 
J- 1 t 



-LJqX 



9AVrm, 



-1 

RMS 1 
''RMS^rVr, 



^3 v RMS^ res 



' l^RMS^ res 

1 - ffAi. 



(18) 
(19) 
(20) 
(21) 



By rewriting the above as two second order equations, we get: 

A 2 = -U) q \2 — T 2 WqAi — gX 3 V RMS V re s 

-T 2 - X A 2 (22) 

^RMS^ res = ^RM 5^rT^es ^ RMS^r^res SX%. (23) 



The above equations can be equivalently found by using: 

Ai - -UqM-u^VxMsgXlVr^-T^Xi (24) 
A 2 = WgAi-T^Aa (25) 



V r , 



= UJ r Z I r es + U r gV RM S^X ~ 7 V res (26) 
= -Z^WrVres- (27) 



Introducing the quantum system voltage and current as: 

V q = VrmsXi, (28) 
h = (-AaTV^flMSWgAs, (29) 

(l24l-(|27]i change into equations for two capacitively coupled 
resonant circuits, with parameters shown in Table |T] 
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